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We discuss the existence and properties of a nontrivial fixed point in /(7?)-gravity, where / is a 
polynomial of order up to six. Within this seven-parameter class of theories, the fixed point has 
three ultraviolet-attractive and four ultraviolet-repulsive directions; this brings further support to 
the hypothesis that gravity is nonperturbatively renormalizabile. 
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The methods of effective Quantum Field Theory 
(QFT) provide an accurate description of low energy 
gravity, including calculable, though currently unmcasur- 
able, quantum corrections On the other hand pertur- 
bative QFT fails when one tries to remove the ultraviolet 
regulator 0. This may signal a failure of QFT to prop- 
erly describe the microscopic features of spacetime, but it 
is also possibile that the problems are only related to the 
use of perturbation theory and that QFT will work when 
more appropriate tools are used. Loop Quantum Grav- 
ity 3 is a nonperturbative approach based on canonical 
quantization methods. Regge calculus [3| and dynam- 
ical triangulations [1] provide discrete nonperturbative 
approximations. Here we pursue a different approach 
that is based on more conventional, continuum, covariant 
quantum field theory methods, and is known as "asymp- 
totic safety" 0. Loosely speaking, a QFT is said to be 
asymptotically safe if there exists a finite dimensional 
space of action functionals (called the ultraviolet critical 
surface) which in the continuum limit are attracted to- 
wards a Fixed Point (FP) of the Renormalization Group 
(RG) flow. For example, a free theory has vanishing beta 
functions, so it is a FP called the GauBian FP. Perturba- 
tion theory describes a neighbourhood of this point. In 
a perturbativcly renormalizable and asymptotically free 
QFT such as QCD, the UV critical surface is parame- 
terized by the couplings that have positive or zero mass 
dimension. Such couplings are called "renormalizable" 
or "relevant" . Asymptotic safety is a generalization of 
this behaviour outside the perturbative domain. 

In the last ten years, evidence for the asymptotic safety 
of gravity has come mainly from the use of the Exact 
RG Equation (ERGE), describing the dependence of a 
coarse-grained effective action functional rfc((f>) on a mo- 
mentum scale k 0, i, i, M, [HI, E [H, H ■ For a review 
see [l^ . This equation has the general form 



+ Rk dtRk 



(1) 



where t ~ log(fc/fco), $ are all the fields present in the 
theory, STr is a generalized functional trace including 
a minus sign for fermionic variables and a factor 2 for 
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complex variables, and Rk is a cutoff that suppresses the 
contribution to the trace of fluctuations with momenta 
below k 0. 

In the spirit of effective QFT, one assumes that rfc[$] 
has the general form 



(2) 



where d [$] are operators constructed with the fields and 
their derivatives that have the required symmetries and 
gi are running couplings of dimension di. Then 



(3) 



where /3j — dtgi- In general the functional ([T]) will con- 
tain infinitely many terms and infinitely many couplings; 
the easiest way of extracting nonperturbative informa- 
tion from the ERGE is to retain only a finite number 
of terms, introduce them in ([1]), evaluate the trace and 
read off the beta functions (3i. Because dtRk in (1) goes 
rapidly to zero for momenta greater than /c, it is not nec- 
essary to use an ultraviolet regulator in this calculation; 
the beta functions are automatically finite. We denote 
gi = gik~'^^ the couplings measured in units of k, and 



A = dm = -digi + Ptk 



(4) 
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where the first term comes from the classical, canonical 
dimension. A FP is defined by the condition f3i = 0. 

The Einstcin-Hilbert truncation consists in retaining 
only terms with up to two derivatives, namely the two- 
parameter Lagrangian density go + giR = Z{2A — R), 
where Z = „ , A is the cosmological constant and G is 
Newton's constant. The existence of a FP in this trunca- 
tion was shown in Q and its properties were extensively 
studied in Q . The stability of the results under changes 
of gauge and cutoff function was interpreted as evidence 
that this FP is not an artifact of the truncation. 

It can be shown in general that if a gravitational FP 
exists, Newton's constant will have an anomalous dimen- 
sion equal to two p^ . in such a way that the propa- 
gator of the theory has effectively a dependence on 
momentum in the UV limit. This suggests that terms 
quadratic in curvature must play an important role. The 
first calculation including such terms was based on the 
three-parameter truncation go + giR + g2R^ [IH. It was 
found that the FP is only slightly shifted relative to the 
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Einstcin-Hilbcrt truncation, and that it is UV-attractive 
in all three directions. Subsequently, the analysis has 
been extended to include also the Ricci squared and Rie- 
mann squared terms, but technical complexities have lim- 
ited the analysis to a one-loop approximation |13|. I n ac- 
cordance with earlier perturbative calculations |l8l. fiot , 
the dimcnsionlcss couplings in the four-derivative sector 
were found to be asymptotically free, while at the FP 
A* = 0.221 and G* = 1.389. The nontrivial FP is attrac- 
tive in all five parameters. The results of [ll| suggest 
that in a more accurate treatment the FP-value of the 
four-derivative couplings would be shifted to a finite and 
nonzero value. 

In addition to the existence of a FP, nonpcrturbative 
rcnormalizability requires that the UV critical surface, 
defined as the locus of points whose trajectories arc at- 
tracted towards the FP when t — > oo, has finite dimen- 
sionality. If this condition is met, the requirement of 
being attracted to the FP, which guarantees a sensible 
UV behaviour, fixes all couplings up to a finite number 
of free parameters that have to be determined by exper- 
iment. This ensures that the theory will be predictive. 

The attractivity properties of a FP are determined by 
the signs of the critical exponents 'di, defined to be minus 
the eigenvalues of the linearized flow matrix 



m 



(5) 



The couplings corresponding to negative eigenvalues 
(positive critical exponent) are called relevant and 
parametrize the UV critical surface; they are attracted 
towards the FP in the UV and can have arbitrary val- 
ues. The ones that correspond to positive eigenvalues 
(negative critical exponents) are called irrelevant; they 
are repelled by the FP and must be set to zero. One can 
show from (4) that at the Gaussian FP -di ^ di, so the 
relevant couplings are the ones that are power-counting 
renormalizable (or marginally rcnormalizable). In a local 
theory they are usually finite in number. 

At a nontrivial FP the canonical dimensions receive 
loop corrections. However, such corrections are expected 
to be finite, in which case at most finitely many critical 
exponents could have different sign from the canonical 
dimension di. Therefore, it is gencrically expected that 
at any FP in a local theory there will only be a finite 
number of relevant couplings Q. 

It is clearly important to substantiate this general ex- 
pectation with explicit calculations. Some evidence for 
the finite dimensionality of the UV critical surface comes 
from the leading order in the l/N approximation, where 
the contribution of graviton loops to the beta functions 
is neglected relative to the contribution of matter fields 
p^ . In this approximation the couplings with three or 
more powers of curvature are irrelevant and the critical 
surface is five-dimensional. In this letter we present for 
the first time evidence that in pure gravity the UV criti- 
cal surface is finite dimensional. 



We will work with the so-called "/(i?)-gravity" theo- 
ries, where the operators in (2) are restricted to be powers 
of the Ricci scalar: Oi = / d'^Xy/gR'^. These theories have 
attracted much attention recently in cosmological appli- 
cations [20| . The quantization of such theories at one 
loop has been discussed in [2l| . Here we analyze the RG 
flow of this type of theories, assuming that / is a polyno- 
mial of order n < 6. In order to apply (1) to gravity we 
choose a background gauge condition; the metric is split 
into a background field and a (not necessarily small) 
quantum field g^^, = g^i, + /i^jy. The (Euclidean) 

action is approximated by 



Tfel^] = Yl / + Sgf + Sc , 



(6) 



where $ = {h^^^c^,Cu\ and the last two terms corre- 
spond to the gauge fixing and the ghost sector 0, fl7| . 
The gauge fixing will have the general form 



Sgf = ^ J d'^xy/gxf,g'"''x^ 



(7) 



where Xu = ^^h^^, — uh^^ (all covariant derivatives 
are with respect to the background metric) . 

The ghost action contains the Fadeev-Popov term 



Sc = j d^x^c^iAghY^c'^. 



(8) 



In the truncation ^ the beta functions can be obtained 
from a calculation of the trace in the r.h.s. of ^ on a 
spherical (Euclidean de Sitter) background. 

As in [3, [m, [l3| we partly diagonalize the propagator 
by using the decomposition 



/iuJ+Vpe^+V^.CM+VMV^cr+-g^^(/i-V2cr). (9) 



The inverse propagator, including the Jacobians due to 
the change of variables ([5]), is given explicitly in [2]| . 
In fill. Il7l| the Jacobians were removed by a nonlocal 
redefinition of the variables ^ and a. This works well up 
to n = 2, but causes divergences for higher truncations. 
For this reason in our calculation we have not performed 
this redefinition. When both procedures are viable, we 
have checked that they lead only to very small numerical 
differences in the final results. The Jacobians can be 
formally exponentiated introducing appropriate auxiliary 
fields and a cutoff is introduced on these variables, too. 

The cutoff operators are chosen so that the modified 
inverse propagator is identical to the inverse propagator 
except for the replacement of z = — by Pk{z) = z + 
Rk{z)] we use exclusively the optimized cutoff functions 
Ri,{z) = (fc2 - z)0{k'^ - z) [13]. This has the advantage 
that knowledge of the heat kernel coefficients up to Bs 
(which contain at most R^ and which we take from [23| l 
is sufficient to calculate all the beta functions. 

At this point, the r.h.s. of the ERGE contains con- 
tributions from the tranverse traceless tensor h'^"^ , the 
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0.00329 
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0.00151 










2 


1.38 


-2.32 


26.9 
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0.00518 


-0.0196 


0.00070 


-0.0097 
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2.71 


-2.27 


2.07 


-4.23 










4 


0.00506 


-0.0206 


0.00027 


-0.0110 


-0.00865 






4 


2.86 


-2.45 


1.55 


-3.91 


-5.22 








5 


0.00507 


-0.0205 


0.00027 


-0.0097 


-0.00803 


-0.00335 




5 


2.53 


-2.69 


1.78 


-4.36 


-3.76-f4.8 


8 i 


-3.76-f4.88 i 




6 


0.00505 


-0.0208 


0.00014 


-0.0102 


-0.00957 


-0.00359 


0.00246 


6 


2.41 


-2.42 


1.50 


-4.11 


-4.42-f5.9 


8 i 


-4.42+5.98 i 


-8.58 



TABLE I: Position of the FP for increasing order n of the 
truncation. 



transverse vector ^, the two scalar components h and 
cr, the transverse and the longitudinal components of the 
ghosts, and the Jacobians. The vector modes correspond- 
ing to the lowest eigenvalue and the scalar modes corre- 
sponding to the two lowest eigenvalues give vanishing h^^ 
in ([9]) and therefore are omitted from the traces. A con- 
siderable simplification comes from choosing the gauge 
p = 0, a = 0: in this gauge the contributions from ^ 
and (J cancel the contributions from the ghosts, leaving 
only four terms, coming from /i^Zj ^ ^^^^ ^he Jacobians. 
Details will be given elsewhere |24l |. 

The differences between our procedure and the one of 
can then be summarized as follows: (i) we use the 
optimized cutoff instead of the exponential cutoffs; (ii) we 
do not redefine ^ and cr; (iii) we remove one extra mode 
from the ghost spectra; (iv) we use the gauge p = 0, 
a ~ Q instead of p = 1, a = 1/Z. These differences and 
the ensuing simplifications allow us to calculate the r.h.s. 
of ([T]) in de Sitter space exactly: it is a rational function 
of R and the couplings gi. The beta functions can be 
extracted from this function by taking derivatives: 



i\ V OR 



■dtVk 



(10) 



R=0 



where V = J d^v^. This has been done using algebraic 
manipulation software, and the limit n < 6 was set by 
the hardware (a standard single-processor machine). 

We can now state our results. Table I gives the position 
of the nontrivial FP and table II gives the critical expo- 
nents, for truncations ranging from n — I (the Einstein- 
Hilbert truncation) to n = 6. For convenience and for the 
sake of comparison with [ll| we also list in Table III the 
FP~values of the cosmological constant, Newton's con- 
stant and their dimensionless product. 

Some comments are in order. First of all, we see that a 
FP with the desired properties exists for all truncations. 
When a new coupling is added, new unphysical FP's tend 
to appear; this is due to the approximation of / by poly- 
nomials. However, among the FP's it has always been 
possible to find one for which the lower couplings and 
critical exponents have values that are close to those of 
the previous truncation. That FP is then identified as 
the nontrivial FP in the new truncation. 

Looking at the columns of Tables I and II we see that 
in general the properties of the FP arc remarkably stable 



TABLE II: Critical exponents for increasing order n of the 
truncation. The first two critical exponents are a complex 
conjugate pair of the form i?' ± ^"i. 





A. 


G, 


A.G, 


n = 1 


0.130 


0.990 


0.13 


n = 2 


0.129 


1.56 


0.20 


n — 6 


0.120 


0.949 


0.114 



TABLE III: Fixed point values of the cosmological constant 
and Newton's constant for n = 1, 2, 6. 



under improvement of the truncation. In particular the 
projection of the flow in the A-G plane agrees well with 
the case n = 1. This confirms the claims made in Q 
about the validity of the Einstcin-Hilbert truncation. 

The greatest deviations seem to occur in the row n = 2, 
and in the columns 32 and • The value of 52* decreases 
steadily with the truncation. The critical exponent 
appears for the first time in the truncation n = 2 with 
a very large value, but it decreases quickly and seems to 
converge around 1.5. This behaviour may be related to 
the fact that g2 is classically a marginal variable. 

As a further test of stability, in the truncation n = 3 we 
have considered the effect of varying the gauge parameter 
a over a wide range of values, keeping p = 0. As a varies 
between and 32, g^*, with i = 0, 1, 2, 3, vary by about 
5 %, I?' changes by 8%, and -d^ change by about 2%. 

When comparing our results for the case n = 2 with 
those of [m , one has to keep in mind that they generally 
depend on the shape of the cutoff function. A signifi- 
cant quantity with very weak dependence on the cutoff 
function is the dimensionless product AG. The value 
0.12 -;- 0.14 given in [ll[ for AG is very close to the value 
we find in all truncations except n = 2. Our value for 
g2* in the n = 2 truncation has the same sign but is be- 
tween one half and one third of their value, depending 
on the cutoff function. This is another manifestation of 
the relatively unstable behaviour of this variable. The 
value given in [ll| for the critical exponent t9' varies in 
the range 2.2 3.2 depending on the shape of the cutoff, 
and is in good agreement with our results, again with 
the exception of the n = 2 truncation. Finally, in [ll| 
the critical exponent t?2 has stably large values of the 
order of 25 with the compact support cutoffs, but varies 
between 28 and 8 with the exponential cutoffs. The val- 
ues at the high end agree well with our result in the 
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n = 2 truncation. The shape dependence that is observed 
with exponential cutoffs can be tafcen as a warning of the 
truncation-dependence of this quantity. 

The most important new result of our calculation is 
that in all truncations the operators from upwards are 
irrelevant. One can conclude that in this class of trun- 
cations the UV critical surface is three-dimensional. Its 
tangent space at the FP is spanned by the three eigen- 
vectors corresponding to the eigenvalues with negative 
real part. In the parametrization (6), it is the three- 
dimensional subspace in defined by the equation: 

33 = 0.00127+0.190.90 + 0.60751 + 1.265 32 

34 = -0.00646 - 0.732 30 -0.0156 31 + 1.880 32 

35 = -0.0155 - 1.132 30 - 0.846 31 + 0.276 32 

36 = -0.0137- 0.59430 - 0.932 31 - 1.283 32 (H) 

Unfortunately, we cannot yet conclude from this calcula- 
tion that the operators Oi with z > 3 would be irrelevant 
if one considered a more general truncation: the beta 
functions that we compute for 3^ are really mixtures of 
the beta functions for various combinations of powers of 
Riemann or Ricci tensors, which, in de Sitter space, are 
all indistinguishable. However, there is a clear trend for 
the eigenvalues to grow with the power of R. In fact, in 
the best available truncation, the real parts of the critical 
exponents differ from their classical values di by at most 



2.1, and there is no tendency for this difference to grow 
for higher powers of R. This is what one expects to find 
in an asymptotically safe theory. 

With a finite dimensional critical surface, one can make 
definite predictions in quantum gravity. The real world 
must correspond to one of the trajectories that emanate 
from the FP and lie in the critical surface. Thus, at some 
sufficiently large but finite value of k one can choose ar- 
bitrarily three couplings, for example 30, 31, 32 and the 
remaining four are then determined by (|lip . These cou- 
plings could then be used to compute the probabilities 
of physical processes, and the relations pT|) . in principle, 
could be tested by experiments. The linear approxima- 
tion is valid only at very high energies, but it should 
be possible to numerically solve the flow equations and 
study the critical surface further away from the FP. 

Extending the results to higher polynomial f{R) trun- 
cations seems to be only a matter of computing power. 
In view of the results obtained here, we expect that a FP 
with three attractive directions will be maintained. 
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